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ABSTRACT 

SN 2008S and the 2008 NGC 300-OT were explosive transients of stars self-obscured by 
very dense, dusty stellar winds. An explosive transient with an un-observed shock break-out 
luminosity of order 1O 1O L0 is required to render the transients little obscured and visible in the 
optical at their peaks. Such a large break-out luminosity then implies that the progenitor stars 
were cool, red supergiants, most probably ~ 9M & extreme AGB (EAGB) stars. As the shocks 
generated by the explosions propagate outward through the dense wind, they produce a shock 
luminosity in soft X-rays that powers the long-lived luminosity of the transients. Unlike typical 
cases of transients exploding into a surrounding circumstellar medium, the progenitor winds 
in these systems are optically thick to soft X-rays, easily absorb radio emission and rapidly 
reform dust destroyed by the peak luminosity of the transients. As a result, X-rays are absorbed 
by the gas and the energy is ultimately radiated by the reformed dust. Three years post-peak, 
both systems are still significantly more luminous than their progenitor stars, but they are again 
fully shrouded by the re-formed dust and only visible in the mid-IR. The high luminosity and 
heavy obscuration may make it difficult to determine the survival of the progenitor stars for 
~ 10 years. However, our model indicates that SN 2008S, but not the NGC 300-OT, should 
now be a detectable X-ray source. SN 2008S has a higher estimated shock velocity and a lower 
density wind, so the X-rays begin to escape at a much earlier phase. 

Subject headings: stars: evolution - supergiants - supernovae: general 



1. Introduction 



SN 2008S (lArbour & Boleslboosh and the 2008 NGC 300-OT ("optical transient". lMonardll2008h were 
initially classified as unremarkable Type Iln transient s that were likely massive star eruptions rather than 



true supe rnovae dSteele et al 



fruitless Jprieto et al 



2008 



Bond et al 



2008). However, searches for optical progenitors proved 



2008bl . iBerger & Soderbergl |2008|), which was peculiar because the transients were 



little obscured and a massive star should h ave be e n easily visi ble. It became clear that they represented a 
new class of massive star transients when IPrietol (12008a) and IPrieto et al.1 (I2008bl) identified the progeni- 
tors in archival Spitzer data as completely obscured stars with luminosities of order 60000L Q and appar- 
ent photospheric temperatures of order 400 K. In mid-IR color-magnitude diagrams, the stars lie at the 
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extreme li mit of the asymptotic branch (AGB) and t hey are extremely rare, with only 1-10 similar stars 
per galaxy jThompson et allbood llGian et allboid) . iThompson etal] J2OO9I) suggested that the M 85-OT- 
2006, SN 1999bw, and SN 2002bu are also likely members of this class, and lKasliwal et all (1201 11) propose 
PTF10FQS as a member. They are clearly not memb ers of the "luminous red nova" population represented 
by V838 Mon and V4332 Sgr ( IKulkarni et all 120071) . Our working hypothesis, which is by no means uni- 
versally accepted, is that the progenitors are extreme AGB (EAGB) stars with masses of order M* ~ 9M© 
cloaked by extremely dense M ~ lO~ 4 M /year winds. 



While a few studies have focused on " traditional" optical examinations of these transients (IBerger et al. 



2009. ISmithet all 2009. S mith et al. 



addition of near-IR and mid-IR data ( Bond et al. 2009. 
20ldlPrieto et alboidlwesson et alboid 



201 Oj), it is impossible to understand even their energetics without the 

Prieto et al.l2009l . lo"hsawa et al 



Bottic ella et al.112009 
Hoffma n et alhoi 



Priet o et al .201 1|). Not only is the near/mid- 



IR emission significant at peak, but it increasingly dominates the emission as time passes, to the point where 
both transients are again invisib le in the optical and near-IR but remain significantly mo re luminous than 
their progenitors in the mid-IR (|Prieto et alJl201CL iHoffman et al.ll201 lL IPrieto et al.ll201 ll) . Thus, one goal 
of the present effort is to provide a self-consistent model of the evolving spectral energy distributions (SED) 
of these transients. 

Unfortunately, there are two possible masses that can be associated with evolved stars having the lumi- 
nosities of the progenitors. They can be cold (3000 K) EAGB stars with M* ~ 9M Q or hotter (> 7500 K) 
~ 15M^ stars, but we hav e no direct information on the progenitor temperatures because of the obscuration. 



In 



Thompson et al.1 (1200 9) we argued for the lower mass scale because extreme AG B stars are expe cted to be 

't ( e.g. Ieo 
all boiOn ■ 



heavily obscured by dusty winds while the more massive stars are not (e.g. iPoelarends et al.ll2008f) . This is 
supported by t he lack of silicate features in the SEDs (e.g. Iwesson et allboid) or in the IRS spectrum of the 



NGC 300-OT (IPrieto et al.ll2009r). since EAGB stars generally hay e graphitic dusts while massive star s have 



silicate dusts (e.g. iGroenewegen et al. I l2009h . 



Umana et al. 



(120101) . further cited bv lSmith et all (I2010h . note 



the presence of PAH-related features in mid-IR spectra of a candidate Luminou s Blue Variable (LB V), and 
thus argue that some more massive stars may also have graphitic dusts. However. IPrieto et al.l (120091) already 
noted that some massive stars have graphitic dust features, but none unite the key points about the NGC 300- 
OT spectrum: the absence of silicate features, the presence of graphitic features and that the wavelengths of 
the graphitic features are indicative of predominantly aliphatic hydrocarbons rather than PAHs, which im- 
plies formation in an environment with little UV radiation. This favors the EAGB model independent of the 
composition. The highest mass stars observed near the transients can be used to estimate an upper, but not 
a lower, bound on the masses of the stars, and there are stars as massive as 15-20M Q near the NGC 300-OT 
(iGogarten et al. 1120091) but not near SN 2008S dPrieto et allboi lh . Because this approach supplies only upper 
bounds, these results again favor the lower mass E AGB stars if the two transients have a common physica l 
origin. Proponents of the higher mass progenitors (IBond et al. 2009, S mith et alj|2009i . IBerger et al.l l2009). 
also favor models in which the progenito r is obscured by a shell of previously ejected material, similar to 
IRC+10420 (e.g. Humphreys et alJll997l). rather than a steady state wind. The lack of mid-IR variability 
from the progenitors (|Prieto et a l. 2008b. IThompson et al.ll2009i) largely ruled out absorption in an expand- 
ing shell since the expansion should have led to an observable drop in the dust temperature over the duration 



-3- 



of the archival observations of the progenitors, again favoring the dense winds observed around EAGB stars. 
We will find that the nature of the transient both constrains the radius of the progenitor star to be very large 
and requires the existence of a very dense wind, further reinforcing the EAGB hypothesis. 



The nature of the transients is also debated. In 



Thompson et al. ( 20091) we outlined a broad range of 



possibilities including stellar eruptions, sub-luminous SN and the final phases of stripping stars to leave 
a massive white dwarf. Potentially the most interesting of these possibilities is an el ectron capture SN 
(ecSN), since they are predicted to be s ub-l uminous and associ ated with EAGB stars (e.g. iPumo et al. Il2009h . 



Bond et al 



(2009) 



Botticella et al. 



(2009) favored the ecSN scenario and 



Berger et all (120091) and lSmith et al.l (12009b firmly fav or non-explosive, stellar eruptions, 



Kashi et al. 



(1201 01) suggested binary interactions. At 



present, the only evidence against explosive transients is the failure to d etect either X-ray or radio emission 
(|Chandra & Soderbergil2008L iBotticella et al.l I2QQ9L iBerger et all 12009). The present day obscuration and 
luminosities mean that the simplest test of these hypothesis, the survival of the progenitor, cannot be carried 
out. While we will be unable to answer this question precisely, we will find the transients had to be explosive 
in nature. 

Our goal here is to combine the available observations of the progenitors and the transients to produce 
a self-consistent model that addresses the puzzles about these events outlined above. We make only three 
hypotheses. First, the progenitors were obscured by dense dusty winds, whose properties we constrain in 
^2] Second, that the transient ejected a significant am ount of mass at veloc i ties characterized by the 500 



1000 km/s velocities ob s erved near the tran sient peak (|Botticella et al.ll2009L iBond et al.ll2009L ISmith et al. 



2009 IBerger et al.ll2009L ISmith et al.l I2010T) . Third, that the dust destroyed by the transient subsequently 
reforms. Section $3] combines these elements to build quantitative models of the transient light curves in 
terms of the evolution of the luminosity, shock radius, characteristic dust radius, and dust optical depth. In 
§@]we examine the physical consequences of the model. In 94. II we show how the destruction of dust needed 
to render the transients optically visible at peak requires an explosive transient from a cool, red supergiant. 
In §14.2l we discuss recombination, cooling and dust re-formation in the wind exterior to the expanding shock. 
In S4.3l we propose that the late time luminosity is powered by X-ray emission from shock heating the dense 
wind that is subsequently absorbed and finally escapes in the mid-IR. We also predict the X-ray emissions, 
finding that while neither source should have been detected in existing X-ray observations, SN 2008S should 
now be a detectable X-ray source while the NGC 300-OT should not. Similarly, we expect early phase radio 
emission to be absorbed. Finally, in 94.4l we discuss the expected evolution of the optical, near-IR and mid- 
IR emissions and in 34. 5 1 the possibilities for detecting surviving progenitor stars. In $5] we summarize the 
results, relate them to other similar transients and discuss future tests of the hypothesis. We adopt distances 
of 1.88 Mpc for NGC 300 (bieren et alJl2005h and 5.6 Mpc for SN 2008S i n NGC 6946 dsahu et alJ boo3) 
and Galactic foreground extinctions of E(B-V) = 0.013 and 0.342 based on lSchlegel et all (11998). 
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2. The Progenitors 



An EAGB star is surrounded by a radiatively accelerate d wind of mass loss rate M and asymptotic 
velocity v w (e.g. IPoelarends et al.ll2008l . Ilvezic & Elitzurl l2010). We are only interested in the radial regime 
where the wind has reached its asymptotic velocity and the density distribution is p gas = M / 'Ai\:v w r 2 . Thus, 
one fundamental parameter is the wind density parameter, M/v w , which is proportional to both the wind 
optical depth at any wavelength and the luminosity produced by an outgoing shock wave (see 94.31) . We fit 
the observed spectral energy di stributions (SEP) of the SN 2008S and N GC 300-OT progenitors using the 
dusty wind models of DUSTY (Ilvezic & Eh^zurl 1 19971 Ilvezic et al.lll999l) . Fig. Q] shows the two progenitor 
SEDs fit as very cool stars (T* = 2500 K) surrounded by dense winds with inner edge dust temperatures of 
1500 K. 

The models are not unique in terms of stellar temperatures or inner edge dust temperatures, but the 
stellar luminosity L* and the density parameter of the wind M/v w are well-defined, as shown in Fig. [2J The 
SN 2008S progenitor had luminosity L* = 1O 4 - 61±O O4 L and requires a wind with log(M/v H ) = -4.93 ±0.12, 
where logM ~ -3.8 ±0.1 is inM Q /year and v w ~ 15 ±4 is in km/s. The mass loss and velocity estimates are 
based on DUSTY's self-consistent wind acceleration models. The NGC 300-OT progenitor had luminosity 
L* = 10 4 - 90±0 04 L Q and requires a wind with log(M/v w ) = -4.26 ±0.06, where logM ~ -3.2 ±0.1 is in 
M Q /year and v w ~ 12 ± 3 is in km/s. The time required for the wind to expand to these radii i s consistent 
with the <^ 10 4 year limit on the lifetimes of the progenitors in this state ([Thompson et al.ll2009l) . 



As emphasized by Ilvezic & Elitzurl (120101) . there is a limit to the mass loss rate that can be driven by 
radiation pressure on dust, which they refer to as the "reddening bound." If these are extreme AGB stars 
with masses of order M* ~ 9M & , the progenitors are near but below this limit for graphitic dust, but above 
the limit for silicate dusts. This adds further support for graphitic over silicate dusts in these systems (see 
Prieto et at]l2009 . Wesson et alJboiol). The winds of the se systems are still extreme, and not those of typical 
AGB stars with dusty winds (e.g. iMatsuura et al.ll2009l) . But we a lready know that these are very atypica l 



stars, with short lifetimes and only a few such systems per galaxy ( Thompson et al.ll2009llKhan et al. 



The parameter which best scales our later results is the radius at which the visual optical depth of the 
progenitor wind is unity, R(tv = 1). The DUSTY model estimates of this radius are shown in Fig. and 
they can be well-approximated by 

log/?(ry = l) = 16.92 + 2.24(logL,-4.61)-a(logL,-4.61) 2 (1) 
log/?(r y = l) = 17.57 + 1.32(logL* -4.90) 

for SN 2008S and the NGC 300-OT, respectively. Here, a = 16 if logL* < 4.61 and a = if logL* > 4.61, 
and the luminosity ranges are 4.51 < logL* < 4.79 (SN 2008S) and 4.78 < logL* < 5.03 (NGC 300-OT). 
This radius is directly related to the wind density parameter 

M _ 4ttR(t v = 1) 

with a visual opacity Ky = 137 cm 2 /g, as shown in Fig. [2] The wind d ensity depends on bo th the bulk density 
of the dust, where we have used pi, u ^ = 2.2 g/cm 3 for graphitic dust ( Ivezic & Elitzur 201ol) . and the gas-to- 



(2) 
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dust ratio, where we are presenting results for DUSTY's default ratio of r gc j = 200. The density parameter, 
M/v w , mass loss rate M and wind velocity v w scale as r g dPbuik, (fgdPbuik) 3 ^ 4 and (r gc ipbuik)~ 1 ^ , respectively. 
More importantly, the models of the progenitor's SED directly constrain only R(ry = 1) oc Mny /v w , so we 
can increase the gas density of the wind M/v w oc r gc i by reducing the dust opacity per unit mass Ky oc r~j, 
while leaving the dust radiative transfer unchanged and ignoring the problem of accelerating the wind. The 
default DUSTY value of r ?l -i = 200 is larger than the ratio typical of the ISM (r g( j ~ 100), and this is generally 
true of AGB star winds (|Ivezic & Elitzurll2010i) . While we leave the gas-to-dust ratio fixed to this default, 
it does provide a means of changing the gas density and thus the shock luminosity and X-ray emission we 
discuss in ^43] without altering the models for the SEDs. 



Most models of these systems have tried to normalize t he parameters using a dust photosphere at 
radius Rh* with temperature Th*, where = AttR^ctT^ (e .g. IPrieto et al.ll200 8b. Tho mpson et al.1 12009 . 



Wesson et al. 



20101) fit dusty radiation 



Botticellaet all 120091. IPrieto et all 120091. iBerger et all 120091. while 
transfer models). There is no well-defined mid-IR photosphere, however, because the optical depth varies 
rapidly with wavelength. In the DUSTY models, the ratios between the visual opacity and those at 5, 7 
and 10/im are Ky / k$ ~ 79, Ky / k-j ~ 138 and Ky / k\o ~ 200, respectively. If we simply estimate where the 
radius R(tj = 1) at which the 7jum (chosen to lie close to the peaks of the two SEDs) optical depth is unity, 
we can illustrate the simple scaling between radius and wavelength-dependent optical depths. The simple 
p oc 1/r 2 model prediction that R(t 7 ) = R(t v )(t v /t 7 )(k 7 / k v ) = 10 14 78 and 10 15 42 cm for r 7 = t v = 1 and 
kv/kj = 138.5 agrees well with the DUSTY models, as illustrated in Fig. [3] 



3. The Evolution of the Transients 



The next step is to characterize the evolution of the transients. Figs. 0] and [5] show snapshots of the 
evolution of the SEDs of SN 2008S and the NGC 300-OT, and Figs.[6]an d|7]show the evo l ution o f the optical, 
near-I R and mid-IR l i ght curves. The data for SN 2008S comes fromlBotticella et all (120091) . Ismith et al. 



(2009), 



Wesson et al. 



(2010), 



comes fromlBond et al 



Prieto et al 



(2010) and Hoffman 



et al.ll 



(200% Berger et all (12009). 



Prieto et al 



4201l|)_and Hoffma n et all (1201 II). while that for th e NGC 300-OT 



2011 



(120091) . IPrieto et all (120101) . lOhsawa et al. 



). We will model these SEDs with DUSTY, ignoring light trave l times ("dust 



echoes", Wright 1980, Dwek 1983), as these have already been considered in some detail by IWesson et al. 
(|2010i ). and we are primarily interested in the later time behavior when they are unimportant. This does 
mean that our very early time models will somewhat mis-characterize the dust properties, where the general 
sense is to underestimate the radius of the dust and its optical depth. However, on the logarithmic scales in 
which we are interested, these are minor corrections, and we will incorporate echoes into our light curve 
models to explore their effects. We can see from Figs.|4]and[5]that the transients, like the progenitors, cannot 
be well-characterized without mid-IR observations. In fact, in the most plausible models, the NGC 300-OT 
peaked in the mid-IR rather than the optical. In both cases, as the lum inosity drops, the SEP shifts from th e 
optical to the mid-IR and the transients again become fully obscured (IPrieto et al.ll201(X IPrieto et al.ll201ll) . 



The relevant variables for describing the transients are the (bolometric) luminosity and temperature of 
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the transient and the characteristic radius and optical depth of any surrounding dust. We first estimated these 
quantities using DUSTY models of the individual epochs shown in Figs.[4]and[5] DUSTY's structure makes 
it difficult to produce an evolving series of models to fit the light curves, so we fit the light curves using sim- 
pler physical models constrained by both the DUSTY models and the light curves. In these models for the 
evolution of the light curves, the transients are modeled with a varying luminosity, L(t) and temperature T(t) 
black body transient, obscured by a time varying dust optical depth ry(?) located at radius r^if). The shock 
veloc ities were constrai ned by the estimates of v s = 1100 and 560 km/s for SN 2008S and the NGC 300-OT 



from [Smith et al.1 (I2010f) . The luminosity evolution can almost be modeled as an exponential decline plus 
a constant, L = L exp(-?/? )+^i, and the temperature T(t) can almost be held constant, but this proved to 
be inadequate for detailed models of the light curves. Thus, we modeled the logarithms of the luminosity, 
temperature and optical depths as piecewise linear functions. We also included an initial unobserved lumi- 
nosity spike that heats the dust at radius r^{t = 0) to a destruction temperature of Tdest = 1500 K. The dust 
emission is then time averaged to include the effects of "dust echoes" in the light curve models. The full 
details of the model are presented in the Appendix. Obviously, this model has some ad hoc features, but it 
produces remarkably good fits to the light curves given its simplicity, as shown in Figs. [6] and [7] The modest 
mismatches of the DUSTY and semi-analytic estimates of the dust radii and optical depths illustrate some 
of the biases created by the approximations. 

Fig. [8] summarizes the luminosity evolution of the transients. The luminosities drop exponentially in 
the early phases and then settle at a nearly constant plateau, given the uncertainties in reconstructing the 
luminosity with variable wavelength coverage. If we approximate the luminosity as L = Loexp(-?/?o)+£i, 
we find L ~ 1O 7 ' 3 L , L\ ~ 1O 5 8 L and t ~ 48 days for SN 2008S and L ~ 1O 78 L , L x ~ 1O 55 L 
and to ~ 39 days for the NGC 300-OT. The uncertainties are almost entirely systematic, so we do not 
report formal statistical uncertainties. Thus, the NGC 300-OT transient was actually brighter, rather than 
fainter than SN 2008S, and radiated more energy in the transient phase (£~ 3 x 10 47 ergs for SN 2008S 
versus E ~ 8 x 10 ergs for the NGC 300-OT). Here we have calculated only the contribution from the 
exponential term. The energy associated with the constant term is unimportant on time scales of to, but 
matches it if the emission continues for 4 and 21 years for SN 2008S and the NGC 300-OT, respectively. 
SN 2008S has nearly reached this point, while the time scale is much longer for the NGC 300-OT because 
of its stronger transient and lower plateau luminosity. The inferred peak lumin osities, Lq g 8 x 10 40 and 
2 x 10 41 ergs/sec are comparable to those of low luminosity Type IIP SNe (e.g. IPastorello et all 120041) . In 
the plateau phase, SN 2008S is brighter than the NGC 300-OT, and both transients remain significantly 
brighter than their progenitor stars. In both cases, the energetics and luminosities derived from traditional 
optical approaches are highly misleading - for example, the luminosity history of these transients including 
the mid-IR emissi ons is very different from the histories based only on optical observations presented by 



Smith et all (12010). 



Fig. [9] shows the evolution of the optical depths. Where the progenitors had enormous dust optical 
depth s, the transients at peak had little obscura t ion, only Ty — 0.3 for SN 2008S and ry ~ 2 for the NGC 300- 
OT (|Prieto et al.ll2008bL botticella et allbood bond et allbood Wesson et allboiOj) . Note the congruence 



between the optical depth estimates at peak in Fig.|9]and the optical depth of the progenitor wind outside the 
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estimate of the initial radius of the dust in Fig. [TOj As the transients evolve, the dust optical depth steadily 
rises, with log Ty ~ 1 .9 ± 0.7 for the final DUSTY epoch in both cases, where the large uncertainties are due 
to the incomplete wavelength coverage in the final epochs. The optical depth does not, however, initially rise 
to the total optical depth of the progenitor wind outside the estimated emission radius, but it does appear to 
steadily approach such optical depths at later times. As the dust destruction radius recedes, the dust appears 
to reform from the outside inwards, but some time is required before the optical depth returns to the full pre- 
transient values. The time scale for restoring most of the dust opacity outside the expanding shock appears 
to be approximately 2-3 years in both systems. 

In the DUSTY models, the preferred transient radiation temperatures drop from 7000-8000 K in the 
first month to 5500-6500 K after 100 days, similar to the evolution in the semi-analytic models. These 
tend to be somewhat warmer than previous models where the optical temperature was estimated by fitting 
the emission without a self-consistent model for the dust absorption. At later times, the heavy obscuration 
makes it impossible to estimate the intrinsic temperature of the transient. As the dust reforms and the 
characteristic radius of the dust shrinks, the characteristic dust temperature rises with a peak temperature 
after roughly 8-9 months. This supports the strong near-IR emission in both sources even as the optical 
emission is collapsing. Fig. \\0\ shows that the initial transient must destroy the dust out to a radius of 
order 10 17 cm, and it then seems to reform from the outside inwards as the dust destruction radius retreats. 
Eventually the radius of the dust approaches the shock radius and the SED models are consistent with the 
dust radius starting to expand outwards again. 

While there is no amb iguity about the pres ence of dust near the peak of SN 2008 S due to the early- 



phase mid-IR detections of IWesson et al.l (|2010h . our c laim of a significan t dust optical depth at the peak 
of the NGC 300-OT differs strongly from the view of iBerger et all (12009b that th ere was no dust o ptical 
depth near the p eak and perhaps we akly with the modest obscuration proposed by iBond et all (120091) . The 
first argument in lBerger et all ((2009) is that their optical/UV SED at 43 days is well-fit by a ~ 4700 K black 
body with no extinction based on the photometry from z-band to the SWIFT UVW1 (2500A) ba nd. In Fig. [5] 
we have extended this wavelength baseline to include the near-IR data from IBond et al.l (120091) and shifted 
to a slightly earlier date (25 May i nstead of 31 May) in order to include the shorter wavelength SWIFT 
UVW2 (1900A) measurement from lBerger et al. ( 2009 ) with less of a correction based on the evolution of 
the UVW1 magnitudes. Compared to a cool, 4700 K black body, the SED has both a near-IR and UVW2 
excess, which the DUSTY models resolve by using a hotter (~ 9000 K) spectrum with an optical depth of 
ry ~ 2.3 ± 0.3. This o pt ical depth is also the amount required to reproduce the Balmer decrement observed 



bv bereer et aD J2009r i. liond et all J2009h interpret the spectra as suggesting that the photosphere has the 
intrinsic color/temperature of an F star, which then suggests ry ~ 1, but they did not directly fit the SEDs. 



Berger et all (|2009h argue against any surviving dust based on the absence of significant Na I D ab- 



sorption in their spectra. The lack of absorption is remarkable because the total sodium column density 
of the circumsteilar medium is of order 10 15_16 ry cm -2 , where ry <; 10 2 is the visual optical depth of the 
material obscuring the progenitor. Unsaturated Na I D absorption lines mean that all the sodium exterior 
to the expanding shock has to be photoionized independent of the survival of any dust. The flux required 
to photoionize a p oc 1/r 2 wind from r,„ to r out scales as Q = <2o(l ~ r in/ r out) where for pure hydrogen 
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Q n = M 2 an/47 rvl,m 2 ri n , which is very different from the Q oc ij mt scaling of a uniform density medium (e.g. 



Franssonlll982h . In a wind, once there are enough ionizing photons, Q ~ Qo, to ionize the inner parts of 



the wind, virtually no additional ionizing photons are required to ionize the entire wind. For example, at 
100 days in the NGC 300-OT, the shock radius of ~ 10 14 6 cm is well inside the progenitor's dust "photo- 
sphere" at ~ 10 15 - 6 cm, so the sodium in the dense parts of the wind has to be photoionized. Expanding the 
sodium Stromgren sphere from the dust "photosphere" to infinity then only requires ~ 10% more ionizing 
photons. Moreover, while the transients cannot keep hydrogen and helium fully ionized due to their low 
radiation temperatures, they are hot enough to produce sodium ionizing photons (see 94.21 ), thereby allow- 
ing hydrogen and helium to recombine while keeping sodium photoionized. This leads to very long sodium 
recombination times due to the lack of electrons, making it sill easier to keep sodium photoionized. Thus, at 
the phases of the spectrosc opic observations, w e expect no Na I D absorption from the wind, and the weak 
Na I D absorption found by lBerger et all (120091) is due to a small amount of interstellar absorption. 



4. A Physical Model for the Transients 



All these features of the evolution of the transient can be explained as a consequence of an explosive 
transient from a red supergiant surrounded by a dense, dusty wind. We assume a wind normalized by the 
prog enitor properties from $2] and shock velocities near the values of v s = 1 100 and 560 km/s proposed 
by ISmith et all J20ld) for SN 2008S and the NGC 300-OT. In 94J] we consider the implications of the 
dust radius and optical depths shortly after the peak of the transient, and find that it requires an explosive 
transient from a red supergiant. In 94.21 we consider recombination, cooling and dust reformation in the 
progenitor wind following this luminosity peak. In 94.3l we explain the roughly constant late time luminosity 
as emission from an expanding shock propagating through the dense, progenitor wind. Eventually, the 
wind column density outside the shock is low enough for the X-rays to be visible, but at early times we 
expect neither detectable X-ray nor radio emission from the shock. Finally, in 94.41 and 94. 5 1 we discuss the 
evolution of the optical, near-IR and mid-IR emissions and the potential detectability of surviving progenitor 
stars. 



The X-ray and UV emission produced by an SN shock wave propagating through circumstellar ma- 
terial has been considered extensively (e . g. Chevalier 1982 . Fransson 1982 . E^undqvist & Franssonlll988[ 
Ichugailll992l . Ichevalier & Franssor]ll994[[chugai & Danzigerll 19941 ). Most of these models were develop ed 



for normal supernovae with early and long-lived X-ray emission (e.g. SN 1978K, Schlegel et alJll999l or 
SN 1980K. ICanizares et al.lll982h . Our outline of a model for SN 2008S and the NGC 300-OT are simple 
versions of these earlier models. The resulting expectations only differ from these models because these 
transients have significantly lower luminosities and shock velocities and are expanding into exceptionally 
dense, dusty winds. 
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4.1. Dust Destruction, Shock Breakout and the Progenitor Radius 



In addition to showing the estimates for the typical radius of the dust emission, Fig. [10] also shows the 
radii corresponding to the dust "photosphere" of the progenitor at 7^m, R(tj = 1), and the radius where the 
visual optical depth becomes unity, /?(ry = 1), as compared to the dust destruction radius predicted by the 
simple exponen t ial plu s plateau models for the bolometric lumin osity shown in Fig. B p As already noted by 



Botticella et all (l2009h . IPrieto et all (|2009), 



Wesson et al 



(|2010h and bhsawa et all (|2010h . and illustrated 



by our models in Fig.[lOj the surviving dust lies in the un-shocked wind material and so can only have been 
destroyed radiatively. Moreover, the required dust destruction radii are an order of magnitude larger then the 
destruction radius predicted from the observed transient luminosities, as also shown in Fig. [TO] The observed 
transients could only destroy the dust out to radii near the mid-IR "photospheres" of the progenitors, which 
would have left the transients heavily obscured at peak with T pea k > 10 instead of Tpeak ~ 1. 

This leads us to the conclusion that the transients must have had short, un-observed, high luminosity 
spikes, and were thus explosive transients. If T pea k is the visual extinction at the peak of the transient, then we 
must destroy the dust out to the radius Rj est — T P eak^ T v - 1) where R(t v = 1) is the radius in the progenitor 
wind where the visual optical depth was unity. The luminosity needed to destroy the dust out to radius Rdest 
is 

2 



Lpeak ^ 167rcr T dest R dest Q rat = 4 x 10 1 r v L eak 



' QlltTdeA ( R(ry = 1) 
1500 K ) V 10 17 cm 



(3) 



where T^ est — 1500 K is the dust destruction temperature and Q rat = Q(T c i es t) / Q(T pea k) is the ratio of the 
Planck-averaged absorption effi ciencies at the temperatures T pea k of the break out radiation and Td est (e.g. 
Dwekll985[braine & Leell984n . Q rat depends on grain size because the smaller grains reach higher temper- 
atures for a given radiation field, so for a more accura t e esti mate of the requi red luminosities we d etermined 
the fraction of the dust destroyed for a iMathis et al] (| 19771) distribution of IDraine & Led (|1984h graphitic 
dust with the 0.005 < a < 0.25/im size distribution used by DUSTY. A transient peak with luminosity L peak 
and radiation temperature T pea k = 50000 K destroys 50% of the dust opacity to radius 



Rd50 



-5/4 



(4) 



and 90% of the dust opacity to R^go c± Rdso/2. Since we require Rdso — R(tv = 1) / 1 Tpeak — 10 cm, we must 



have 



^peak,50 



3 x 



\TpeakW Cm J 



(- 



l dest 



2.5 



(5) 



V 1500 K 

with the luminosity L peak ^ needed to destroy 90% of the opacity being roughly 4 times higher. Such a high 
peak luminosity can only be explained by the luminosity produced from a shock breaking out of the stellar 
photosphere, which in turn requires that these were explosive transients rather than some form of radiation 
pressure driven eruption. 

Producing such a high luminosity with a shock breakout also implies that th e progenitor star s were cold, 
red supergiants. We show this following the approximate breakout model of lOfek et all (120101) (following 
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Falk & Arnettlll977l IWaxman et all 120071) . The shock breaks out when the radiative diffusion time scale 
matches the expansion time scale, = i^tPs^/c = t exp = r s /v s , where kj ~ 0.4 cm 2 /g is the Thomson 
opacity. This determines the density p s of the breakout region in terms of the break out radius r s , 

Ps = ^—- (6) 

v s K T r s 

The post-shock region is radiation dominated, so the radiative energy density has to match the energy density 
of the shock, aT^ = (7 /2)p s v 2 s , where we have assumed the density compression for T = 4/3. Thus, 

Off-l^-. (V) 
By setting the break-out luminosity = L pea k,50> we find that the radius is 

, ,„af R(rv=l) \V 71„, \ 5/2 /1000km/s\ 



This is relatively crude estimate, so we have introduced the dimensionless factor / to monitor the effects of 
the approximation on other estimates. For example, / = 0.5 if we require Lj, ea fc,90 instead of L pea ^$Q. If we 
interpret this as the shock breaking out of the wind, then the wind density parameter must be 



M - 47rcr -v ^- 10 -3.2 ( *(7V = 1) V( T Jest y /2 / lOOOkm/s \ 2 M /year 
v w v s k t \TpeakW 1 era) \1500Ky \ fv s ) km/s 

which is an order of magnitude larger than the wind densities inferred from the progenitor properties (Fig. O. 
Thus, the shock breakout radius must correspond to the stellar radius R* ~ r s . This implies a stellar temper- 
ature of 

r. . f-=V) " 4 K 4000 (-£-) " 4 (W)"-') (wo*) 5/4 ( fi. ) " 2 k. (10 ) 

For our nominal estimates of the progenitor luminosities and wind properties from $2] and the optical depths 
at peak from $31 this yields stellar temperatures of order 1400 K with formal uncertainties of order 30%, 
but the estimate is really dominated by systematic uncertainties in the models and dust destruction physics - 
temperature is a variable where logarithmic accuracy is not really adequate. Nonetheless, the estimate is far 
more consistent with a cool, red supergiant as the progenitor, an EAGB star, than with a hotter, bluer star. 
The break out luminosity increases with the size of the star, L s oc r s , because the larger surface area (oc r 2 ) 
matters more than the lower shock temperature (r/ oc 1 /r s , Eqn. |7]). The duration of the break out peak is 
of order the light crossing time r s /c ~ 10 4 sec or less, so the energy released in the peak is less than that of 
the observed transient. 



4.2. Recombination, Cooling and Dust Reformation 



The transient models from $3] appear to show that (1) the dust optical depth rebuilds from the outside 
in, (2) that some optical depth returns quickly, but (3) that it takes several years for the optical depth to 
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rebuild to levels close to that of the progenitor wind. Sadly we lack the continuous, detailed monitoring of 
the mid-IR spectral energy distributions that would have allowed a careful exploration of this process. In 
this section we outline some of the relevant physical issues and Fig. [TT] outlines our model for the process 

The temperature of the shock break out radiation is high, T s ~ 50000 K (Eqn. |7]), so the wind is pho- 
toionized by the break out radiation. For pure hydrogen, the minimum ionizing flux needed to maintain the 
ionization is 

4na a R(T V =lf _ 4 , / a, \ 1 140 cm 2 /g \ 2 / fifty = 1) \ ( R(tv = \) \ 



which cannot be supplied by the transients even at their observed peaks (low luminosity, low temperature, 
and R(ry = l)/R ~ 10 3 ). We scale this and the subsequent expressions to the radius R = 10~ 1 R(t v = 1) ~ 
10 16 cm, roughly corresponding to the shock radius at the present time (~ 3 years at 1000 km/s). Thus, the 
wind recombines on time scales (for hydrogen) of 

= , mp R(r v = l )Kv / /lO-iW/sX /_J0R_y 

(12) 

that are very short in the interior and only become relatively long as you approach R(ry = 1), as shown in 
Fig. [H] While the transients cannot keep hydrogen and helium in the wind fully ionized due to their low 
temperatures, they have no difficulty keeping sodium photoionized, as discussed in $3] At late times, as 
we discuss in §14.31 the source of the luminosity is non-thermal emission from a shock, and ultimately this 
should reionize the wind because the required number of ionizing photons steadily diminishes. 

Then as the gas cools, the dust can reform. The cooling time scale is of order 

3kTum p K V r 2 n nnn ( T \ ( 10" 25 ergs cm 3 /s\ (R{t v = 1)\ / 10/? \ 2 
tcool = P „ * 0.089 -— - f '- years (13) 



2AR(t v = 1) ' V 104 K / V A J \ 10 17 cm J \R(t v = 1) 



where a cooling rate of A ~ 10 25 ergs cm 3 /s is approximately correct for 10 3 < T < 10 4 K (e.g. ISchure et al. 



2009). Thus, the gas also cools rapidly after the peak, with the cooling radius being only slightly smaller 



than the recombination radius in Fig. [Til 

Once the temperature is low enough, T <> Tdest either the dust begins to form anew following a burst 
of nucleation in the now supersaturated vapor or any residual grains can begin to grow. Given a steadily 
declining temperature, nucleation is not a bottleneck, although the radius for inhibiting dust formation is 
larger than the dust destruction radius because less flux is needed to destroy the smallest grains. For a 
7000 K transient and T^ esl = 1500 K, graphitic dust can begin to form outside the radius 

/ r \ 1/2 
,15.9 I L 
\ mRi 

-0 



Rform ~ I0 10 y I I cm, (14) 



estimated from the heating of the smallest grains. We used this estimate of the destruction radius in Fig. [10] 
to illustrate where dust can begin to reform. For both transients, Rf orm is outside the shock radius only for 
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the first ~ 100 days, as shown in Fig. [TT] based on the exponential plus constant luminosity models for the 
transients from $3] For comparison, Fig. [TT] also shows the radius at which 50% of the dust opacity will 
survive aL pea £ = 4 x 1O 1O L break out shock based on Eqn. [5] 

The bottleneck for reforming the dust is the growth rate of the particles after nucleation, where the time 
scale to grow a (spherical) grain of radius a by adding monomers assuming a sticking fraction of unity is of 
order 



^^^^'^^J^Jl-T^jlw^) years (15) 

where v c ~ v w ~ 10 km/s is the collision velocity, and X c is the carbo n mass fraction in t he wind (we have 
used a relatively high value appropriate for a carbon rich EAGB wind (lAbia & Isernll2000l) ). If the particles 
are growing by coagulation rather than monomer accretion, the growth rates are a four times faster, but 
significantly shorter growth times could only come from incorporating significant numbers of the far more 
abundant hydrogen atoms, having a reduced dimensionality (e.g. sheets) or density, or starting from more 
and larger surviving particles at larger radii. In practice, the latter effect is certainly present, so Eqn. Q3] 
does exaggerate the difficulty of reforming dust at larger radii. We should al so note that there is no need to 



fully repopulate the 0.005 < a < 0.25/xm size range of the standard DUSTY iMathis et all (119771) models to 
recover most of the opacity. DUSTY models using distributions with maximum grain sizes of 0.01, 0.025, 
0.05 or 0.1/im instead of 0.25/xm have visual opacities smaller only by factors of 3.4, 3.4, 2.9 and 1.6, 
respectively, consistent with the rapid return of a significant fraction of the pre-transient opacity followed 
by a slower rise towards the total pre-transient opacity. 

Fig.QTJillustrates the cycle of dust formation for SN 2008S. Assuming growth begins after both recom- 
bination is completed and the radiative heating is sufficiently low, dust begins to reform at an intermediate 
radius where the formation time scale is relatively long. However, since the cooling and growth time scales 
are significantly shorter at smaller radii, the optical-depth weighted radius of the dust will rapidly shrink, 
roughly tracking Rf orm , until it reaches the expanding radius of the shock. Our models essentially represent 
the evolution of this optical-depth weighted radius. A variant of the models where we used two shells, one 
at large radius to represent the dust surviving the explosion, and one expanding with the shock, could not 
reproduce the light curves as well. We attempt an explicit model of the evolution of the optical depth in 
Fig. El In this model, the dust begins to reform at each radius after the longer of either the recombina- 
tion time (Eqn. \Y2\ or the time needed to lie outside the region radiatively heated to Tdest (Eqn. [T4l . The 
maximum particle size a max is then determined from Eqn. [I5j and the opacity is found by the scaling of the 
DUSTY model opacities with a max . Qualitatively this reproduces the behavior of the optical depth in the 
models, which are also shown in Fig.Qj] Quantitatively, the optical depth rises too quickly and is somewhat 
low at the time of our last SED model. The early time differences could be due to oversimplifying the 
temperature of the dust at smaller radii where there may be additional heating and destruction mechanisms 
associated with the ionizing radiation from the expanding shock front (see S4.3I ). Overall, however, the 
physical characteristics of the winds are broadly consistent with reforming the necessary amounts of dust. 
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4.3. X-ray Emission from the Shocks Powers the Late-Time Light Curves 



Even after almost three years, these sources are still an order of magnitude more luminous than their 
progenitors, and the present emission is most likely due to shock heating th e dense wind. For a shock wave 
propagating outward at velocity y t » v w , the shock releases energy (e.g. IChevaherll 19821 lFranssonlll982 . 
Chugailll992i Ichusai & Danzigerill994 Ichevalier & Franssonlll994) 



Mv 2itR(t v = l)v; 



2 v v , 



1.2 x 10 6 



Ky 



1000km/s 



3 ( R(T V =l) 

V io 17 



cm 



140 cm 2 /g 



Ky 



L, 



•0 



(16) 



although this energy can be lost due to expansion rather than radiated. In detail, there is both a forward and 
a reverse shock, separated by a contact discontinuity. Models of X-ray emission for typical SN focus on the 
reverse shock because of its lower temperature. Here, the low velocities make the forward shock a source 
of very soft, readily absorbed emission as well. Exactly what should be used for v s is unclear, since the line 
widths in normal Type Iln are not believed to be representative of the true shock v elocity, but a s wit h the 
light curve models in §[3] we will scale results based on the FWHM estimates from I Smith et al. I bold) of 
v s = 1 100 km/s and 560 km/s for SN 2008S and the NGC 300-OT, respectively. Because of the low shock 
velocities, the post-shock temperatures are low 



E s = —m p v 2 = 1.2 



lOOOkm/s 



keV 



(17) 



for a mean molecular weight ^ = 0.6. The optical depth of the swept up wind is the same as that of the 
wind exterior to the shock for the same opacity, so the swept up wind material quickly becomes Thomson 
optically thin, 

R(t v = 1)kt _ nnn ( R(t v =1) \ f 1000 km/s \ / yean 
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(18) 



The low Thomson optical depth of the wind essentially corresponds to the conclusion in 94. II that the wind 
density is too low to be the source of the shock breakout luminosity. Typical models of these shocks consider 
cooling by adiabatic expansion, Compton cooling and free-free emission, but generally consider much faster 
shocks (v s ~ 10 4 km/s) in less dense winds, finding that the early-time cooling is dominated by Compton 
cooling. In these systems, free-free emission dominates over Compton cooling, with a cooling rate of orde r 
A = A n 2 (T /T ) l l 2 for T = 1 keV ~ 10 7 K and A ~ 2 x 10" 24 ergs cm 3 /s (e.g. iRvbicki & Lightman|[l979) . 
so the ratio of the cooling time, t coo i = 3nkT s /2A, to the expansion time, t exp =R/v s ,is 



tcool 



l exp 



0.5 



1000 km/s/ V140 cm2 /g 



Ky 



Ccomp 



10/? 



R(t v = 1) 



(19) 



where C comp > (T+ 1)/(T- 1) = 4 is the ratio of the post-shock density to the pre-shock density. Thus, the 
shock energy will be efficiently radiated as soft X-rays until the shock radius approaches R(ry = 1). 

Unlike typical SN, the combination of high wind densities and a very soft X-ray spectrum mean that the 
X-ray emission is fully absorbed by the unshocke d wind material and rad iated in the mid-IR by the reformed 
dust while the wind is optically thick. Following lDraine & Woods! (119911) . we model the (bound-free) X-ray 
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opacity of the wind as 

k h (E) = 17(£/keV) _3 cm 2 /g (20) 
k z {E) = 85(l + (£/keV) 2 - 5 ) -1 cm 2 /g 

from neutral hydrogen/helium and heavier elements, respectively. This implies X-ray optical depths of 



/140cmVg\ / gfry = l) \ 1 



(E/keV) 2 - 5 ' 

from radius R to infinity. The corrections for the material outside the recombination radius (Eqn. [T2)) are not 
important until very late. The free-free emission spectrum is 

F u Kg(E/E )exp(-E/E Q ) (22) 

where g(x) is the Gaunt factor. For an X-ray temperature of Eq = 1 keV, the fraction of the X-ray energy 
escaping the wind for R/R(tv = 1) = 0.01, 0.1, and 1 is only 0.2%, 5%, and 27%, with mean energies of 
4.8, 2.6 and 1.5 keV. Slightly more escapes for SN 2008S (Eq ~ 1.2 keV) and quite a bit less escapes the 
NGC 300-OT {Eq ~ 0.4 ke V). While soft X -ray absorption from SN shocks has been discussed in the context 



of normal SN models (e.g. iFranssonll 19821 ). the overall spectrum is assumed to be much harder (because v, 



is larger), so only small fractions of the total luminosity are absorbed. 
We can approximate the expected X-ray emission as 

T T 1 — fabs 

Lx — L$ 

1 + texp I tcool 



(23) 



where L$ is the shock luminosity from Eqn. [16] I— fobs i s tne unabsorbed fraction estimated by combining the 
wind model, the X-ray optical depths in Eqn. [22] and the shock energies from Eqn. [IT] and the denominator 
approximates the balance between radiative and expansion losses from Eqn. [19] When the cooling time 
is short, all the energy is radiated, while when it is long, only the fraction t coo i/t exp is radiated and the 
remainder is lost to adiabatic expansion. The evolutions of the expected X-ray luminosity and the mean 
photon energy for the two systems are shown in Fig. [13] Without absorption, both systems should be bright 
X-ray sources. With absorption, however, they are relative faint, hard X-ray sources in th e early phases 



and it is not at all surp rising that they were not detected in the early SWIFT observations dBotticella et al. 



2009, lBergeretal.ll2009h . While SN 2008 S has not been reobserve d, the NGC 300-OT was reobserved 



after roughly two years by both SWIFT dBerger & Chornockl 120101) and Chandra (CXO observation ID 



12238/PI: Williams). Fig. [T3l indicates the count rates at which these various observations would detect 3 
source photons, where we estimated the count rates using the Chandra/SWIFT exposure time calculators and 
absorbed thermal bremsstrahlung models based on the gas column densities in our nominal wind models for 
the two sources. Neither source should have been detectable in any of the observations to date. SN 2008S 
should now be relatively easy to detect, with count rates of order 10~ 3 /sec, but has not been observed, while 
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the NGC 300-OT is essentially undetectable - it is less luminous, softer, and more heavily obscured. Bear in 
mind, however, that at fixed shock luminosity, the detectability of the X-ray emission depends exponentially 
on the shock velocity because the X-ray optical depth of the wind scales as L s /v A s (three powers of v s from 
the shock luminosity, Eqn. [16J and one power because of the slower expansion through the wind), and the 
typical energy scales as v 2 s (Eqn. [FT]). Small reductions in v s dramatically increase the overall absorption. 



Such shocks are also expected to produce radio emission, which has not been observed (iChandra & Soderberg 



2008. iBerger et al.ll2009l) . and this is also used as evidence against an explosive transient. In fact, no early 
time radio emission is likely to escape the wind exterior to the shock because absorbing the X-ray emission 
will create an ionized layer in the unshocked wind exterior to the shock. If this ionized layer has a thickness 
AR, and electron temperature T e ~ 10 4 K, then the free-free optical depth of the layer is 

1.35 
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and these systems will be very optically thick to GHz radio emission in their early phases and should not 
have been detectable. The time scale for the optical depth to fall to t v = 1 is 
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(25) 



ignoring any other potential source of opacity at these wavelengths. Thus, at late times these systems may be 
detectable in the radio. More generally radio emission from such high density winds or from shocked, mas- 
sive shells of material will be very different from the significantly lower re gimes of circumstellar densities 
that are usually considered for radio supernovae (e.g. IChevalier et al. 2006). 



4.4. Mid-IR, near-IR and Optical Emission 

The current mid-IR emission is simply the absorbed X-ray flux combined with any emission from dust 
heated and destroyed by the shock. The absorbed X-ray emission leads to a mid-IR luminosity of 

T T fabs 



(26) 



1 ~^~texp I tcool 

which is simply the absorbed counterpart to Eqn.[23j as long as the dust optical depth remains high. Here we 
are ignoring the extra emission fro m dust directly heated by the shock, which typically radiates ~ 10% of 
the shock luminosity ( Drainefl981 ). The dust temperature is simply determined by the available luminosity, 



since emission by the dust is the only means of radiating the absorbed energy. Thus, the dust temperature 
must be or order 

-fe) 1/4 ^(^)' /4 (") ,/2 (~) 1/2K - - 

where r c j is the dust radius. Driven by the declining luminosity and the increasing dust radius, the mid-IR 
emission slowly declines and shifts to longer wavelengths, as we see in Figs. [6]and|7J This phase lasts until 
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the shock expands to the point that the optical depth of the reformed dust begins to clear. If all the dust has 
reformed, it takes 



t(r x ) = 
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for the optical depth down to the shock to be t\. In Figs. |6] and |7]fhis leads to the sequential recovery of the 
near-IR and then the optical fluxes, with the present period being the phase where the absorption is worst. 
In these models, the near-IR fluxes should begin to recover in the immediate future. 

Unfortunately, the late time details of the light curves in Figs. [6] and |7] are incorrect because the models 
simply continue to assign the transient a black body SED of fixed temperature. This approximation is 
adequate when the optical emission is dominated by the transient rather than the shock or when the emission 
is all absorbed and re-radiated in the mid-IR, but it is not correct when the dust opacities become low enough 
for the line emission from the shock to escape. The SED of the shock will be a complex line spectrum (e.g . 
Fesen&Matonicklll994l for S N 1980K) extendin g from the UV into the mid-IR (e.g. Ullen et alJl2008T) 
rather than a black body. The [Alien et al.1 (12008!) shock models do not reach the wind densities needed 
for these systems, but to the extent we can extrapolate from the highest density models in their survey, 
the dominant emission is in the UV (Lya), followed by the optical (Ha) and mid-IR, with relatively little 
emission in the near-IR. Hence, the extrapolations shown in our light curve models are probably over- 
estimating the recovery in the near-IR. Fig. [14] shows the effect on the light curve of SN 2008S if we model 
the shock emission using the emission l ines and line ratio s of the highest density (n = 10 3 cm -3 ), solar 
metallicity, v s = 1000 km/s models from lAllen et al.l (120081) . The true density is much higher, and so the 
actual line ratios will be strongly affected by the collisional de-excitation effects that limited their survey of 
shock emission spectra to n < 10 3 , but it illustrates the potential differences from the model used in Fig. [6] 
As a result, these shock models have no significant recovery of the near-IR flux as the optical depth drops. 



4.5. Did the Stars Survive? 



Because the shock luminosity is significantly greater than that of the progenitor, we can only indirectly 
constrain the existence of a surviving star. One probe of this ques t ion is the mass of the ejecta. Using the 
expansion velocities of v s = 1 100 and 560 km/s from ISmith et all (120101) for the transients, and assuming 
the radiated energy estimates from $3] are fraction / of the kinetic energy, the ejected masses are of order 
M e ~ O.O2/ _1 M and O.25/ _1 M for SN 2008S and the NGC 300-OT transient. Since destroying the dust 
required an explosive tra nsient, we would expect / ~ 0.01 as is typical of supernovae rather than / ~ 1 
similar to rj Carinae (e.g. ISmith et al.ll2010h . but this argument is too crude to determine the survival of the 
star. The roughly constant shock luminosity means that the ejected mass must significantly exceed the swept 
up wind mass, so M e <; 3Mv s t/v w , which after t ~ 3 years corresponds to M e <; O.O5M and O.25M for 
our standard parameters, weakly requiring / <; 1. Unfortunately, the time scales for estimating M e from the 
evolution of the transients are long. 



This leaves the problem of directly detecting a surviving star. As the shock expands, the pre-existing 
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wind will stop cloaking the progenitor and the shock luminosity will fade due to the decreasing radiative 
efficiency and will likely be dominated by X-ray and line emissions that are not characteristic of the pro- 
genitor. Fig. [14] shows the effect of including the SN 2008S progenitor as a 7* = 3000 K black body in 
the shock model for the late-time emissions. This assumes that no dust has formed in or behind the shock 
but includes the obscuration from the wind exterior to the shock. Obviously, if dust forms either near the 
contact discontinuity between the ejecta and the wind or in the ejecta, there is easily enough material to 
fully cloak the star. Moreover, the luminosity of the shell plus any surviving star is insufficient to maintain 
a temperature above the dust destruction temperature in the region interior to the reverse shock. 



5. Summary 



The available data for the SN 2008S and the NGC 300-OT are consistent with a relatively simple, 
self-consistent physical picture of the systems. The progenitor stars are luminous, cool, red supergiants, 
most likel y EAGB stars based on their positions i n mid-IR CMDs in a short-lived state with very high mass 
loss rates (|Thompson et al.ll2009L iKhan et al.ll2010l) . The very dense winds form so much dust that they are 
optically thick in the mid-IR with dust photospheres on sc ales of a few 10 cm and optical pho tospheres 
near 10 17 cm. They are best modeled with graphitic dusts (IPrieto et al. 2009, Wesson et al.ll2010l) . and the 
radiatively driven dusty wind models of llvezic & Elitzun (120101) can produce such winds at the ~ 9M Q mass 
scale of EAGB stars assuming graphitic, but not silicate, dusts. 

Both stars then underwent an explosive transient with a short duration, very high luminosity shock 
breakout spike. This unobserved luminosity peak is necessary because the observed luminosity peaks are 
not nearly bright enough to destroy dust to the dist ances needed t o hav e the transient peaks little extincted. 
We know from early-time mid-IR observations of Wesson et al. ( 201dh that SN 2008S had a modest sur- 
vivin g dust optica l dep th, and o ur models of the early-time UV through near-IR SED of the NGC 300-OT 



from iBerger et all (120091) and iBond et al.1 (120090 i mply a surviv i ng vis u al optical depth of a few, which is 
consistent with the Balmer decrement observed by iBerger et al.1 (120090 . iBerger et all (120091) argue that the 



weak Na I D absorption in spectra of the NGC 300-OT imply no surviving dust, but we argue that sodium 
must remain photoionized independent of any surviving dust and would be uncorrected with any surviving 
dust if present. Similarly, the dense circumstellar wind is initially opaque to X-rays and requires very lit- 
tle ionization to be opaque to radio emission, so the failure to detect the systems at both X-ray and radio 
wavelengths at early-times is also a natural consequence of the progenitor properties even for an explosive 
transient. The stellar radius required by the shock break out model also requires a large, cool star, again 
consistent with the EAGB picture for the progenitor stars. 

The evolution of the transient SEDs with time is largely controlled by the reformation of the dust in 
the progenitor wind outside the expanding transient shock wave. The average dust radius appears to move 
inwards from near the original optical photosphere as the dust destruction radius shrinks. The dust radius 
then appears to begin to expand again to track the estimated radius of the shock wave passing through the 
wind, although the evidence for this is less striking than the initial shrinking of the dust radius. Initially, the 
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optical depth is significantly less than that of the progenitor wind outside the same radius, but the fraction 
of the dust that reforms steadily rises, with most of the dust having reformed after about 2-3 years. The 
dust radius shrinks faster than the transient luminosity drops, so that the dust temperature peaks 8-9 months 
after the transient peak, and this rising dust temperature balances the falling luminosity to produce relatively 
extended, slowly varying near-IR light curves. The optical fluxes drop very rapidly from the combination of 
the dropping luminosity and the rising optical depths. Although the luminosity eventually becomes roughly 
constant, the optical depth continue s to rise and the systems fade in th e near-IR. At present, neither system 



is detected at optical wavelengths (IPrieto et al 



2010, 



detectable in the near-IR even with HST (IPrieto et al. 



Priet o et al.ll201 11) . and SN 2008S is at best marginally 



201 lb . In the mid-IR, however, both systems are easily 



detected with significantly higher temperatures and luminosities than the progenitor stars. 

The average dust radius in our models is an optical depth weighted mean radius, and the true evolution 
is more complicated. Models with a surviving shell of dust and a shell of forming at the the expanding shock 
radius work poorly - successful models require dust reforming in the exterior wind and are incompatible 
with a shell ejection scenario, just as was true of the pre-transient mid-IR variability. The optical depth 
of the reformed dust is dominated by the smallest radius at which dust can reform because the particle 
growth rates are proportional to density, so the effective dust radius shrinks rapidly, tracking the decaying 
luminosity with a delay, until it encounters the expanding shock. The simple model of a smoothly shrinking 
dust radius elides over a gap at 10 16 to 10 17 cm, between the surviving and reforming dust, where the particle 
growth times are too lon g for rapid reformati on. This pattern o f a receding dust rad ius is also reported for 
the Type IIP SN 2007it (Undrews et alJboi lb and SN 2007od (Undrews et allboiob . although the implied 
optical depths are much lower. Re-forming CSM dust sufficiently quickly requires the high densities of slow, 
high M winds or massive shells of ejecta. In most systems, only the normal scenarios are relevant because 
the particle growth rates are too slow: dust must either survive the explosion, form in cooled shocked CSM 
near the contact discontinuity or form in the un-reverse-shocked ejecta. The high CSM densities that allow 
rapid dust reformation are yet another new (or odd) property of these transients. 

The present day, roughly constant luminosities closely match the luminosity expected from expand- 
ing shocks heating the progenitor winds. The shocks are Thomson optically thin but radiatively efficient, 
leading to high ~ 10 5 ' 5 L Q , soft (0.5-1.0 keV) X-ray luminosities. Soft X-ray opacities are, however, not 
very different from dust opacities, so the dense progenitor winds simply absorb the shock luminosity and 
the energy is ultimately radiated in the mid-IR, leading to the present-day spectral energy distributions. Nei- 
ther source should have been detectable in the X-ray observations of these systems to date, despite having 
intrinsic X-ray luminosities of order 10 39 ergs/s. As the shock moves outward, the X-rays begin to escape 
and should show a spectrum with strong absorption at 0.5-2.0 keV that will be a powerful probe of the 
wind column density. SN 2008S, which has a lower density progenitor wind and a faster shock velocity, 
should now be relatively easy to detect as an X-ray source, while the NGC 300-OT should still be very 
difficult to detect because of its higher density wind and slower shock velocity. Slower shock velocities 
reduce the X-ray detectability partly because the slower shock expansion rate slows the drop in the exterior 
wind column density and partly because the lower shock temperature leads to a softer more easily absorbed 
X-ray spectrum. The detectability does, however, depend exponentially on the assumed shock velocity. The 
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peak X-ray emission is a balance between the diminishing absorption column density and the diminishing 
radiative efficiency of the shock, peaking on time scales of order 10 years for SN 2008S and much later for 
the more slowly evolving, more heavily obscured NGC 300-OT. Absorbing the X-ray emission will produce 
an ionized region in the unshocked wind just outside the shock, and the resulting free-free opacity will also 
block any GHz radio emission produced by the shock for an extended period of time. Thus, the failure to 
detect early-time radio emission is also expected. 

The balance between dust reformation and shock expansion means that the visual optical depth to the 
shock peaked after roughly 2-3 years and should now be dropping linearly with time. The mid-IR luminosity 
should steadily diminish, because there is less absorption, and shift to longer wavelengths, because the larger 
dust radius and diminished luminosity lead to lower dust temperatures. It may be difficult to cleanly separate 
these trends with only 3.6 and 4.5/xm warm Spitzer observations. Near-IR observations will shortly be able 
to penetrate to the shock front, although the visual optical depth to the shock front will not approach unity 
for over a decade. The optical/near-IR emission from the shock region should be a combination of free-free 
emission, line emission driven by X-ray absorption and shock heated dust. In the evolution to date, the 
luminosity was dominated by the shock emission while the dust optical depths are high enough to allow us 
to ignore these details, which will not be true of these later phases. Unfortunately, until the luminosity fades 
or we can penetrate the dust, it is impossible to determine the fate of the progenitor star. 

The author would like to thank J.F. Beacom, J.A. Johnson, M. Pinsonneault, J.-L. Prieto, R. Pogge, 
K. Sellgren, K.Z. Stanek, D. Szczygiel, T.A Thompson and B.E. Wyslouzil for extensive discussions. CSK 
is supported by NSF grant AST-0908816. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration. 



A. Light Curve Fits 

Here we spell out the details of the light curve fitting procedure. We specify the luminosity and optical 
depth on a grid of times t\ with i = 1 ■ ■ • N and t\ = 0. The transient luminosity = L(ti) and temperature 7J(f,-) 
are determined at each grid point, and the value at time t is found by linearly interpolating logL, between 
bracketing times t{ <t < ti+\ or fixed to L^/Tn for t > t^. The visual optical depth Tyj is determined at grid 
points i = 1 • ■ ■ (N- 1), but fixed to the optical depth of the progenitor wind outside the current dust radius 

R(t v = 1) 

T V (r d ) = — — (Al) 

rdit) 

for t > tff, where the dust radius declines exponentially starting from radius ro on time scale to, and then 
begins to expand outwards at the shock velocity v s , 

r d (t) = r expH/%) + v s t . (A2) 

Thus, the last optical depth grid point is tv,n = Ty(^(^v)) an d represents the point where the optical depth 
of the wind outside the shock has fully reformed. The parameters of the model are then the N luminosities 
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L, and temperatures 7}, the N- 1 optical depth points Tyj (i = 1 ■ • -A 7 - 1) and times ?, (i = 2 • • - AO, the shock 
velocity v s , the dust radius parameters tq and to, the ry = 1 radius of the progenitor wind /?(ry = 1), and Q rat , 
which enters into the estimated dust temperature and is a proxy for the more complicated radiation transfer 
of the DUSTY models. 

The emission has four components: a shock break out of luminosity L pea ^, duration r*/c, and temper- 
ature Tpeak = 50000 K; an optical transient of luminosity Lo = L(t)—L^; a shock with luminosity L$ = L^; 
and, finally, the luminosity absorbed and re-radiated by the dust Lp. We model the optical contributions as 
a black bodies, where it is useful to define 



B v (n,T) = \~ n B v {T) 



I 



dv\~ n B v {T) 



(A3) 



a modified Planck function normalized to unit total luminosity, and let /(A) = t\/tv be the optical depth 
ratio in the DUSTY models between the visual optical depth t v and that at wavelength A, t\. The escaping 
radiation from the transient is then 

L u . (t) = L B^(0, T{t))e- Tv{t)f(X) (A4) 
and the fraction of the optical transient luminosity escaping is 

fa,esc = J dvB„(0,T(t))e- Tvit)fW (A5) 

The shock luminosity is assumed to be initially generated as X-rays at the energy of Eqn. [T7J and then 
absorbed by the wind exterior to the shock using the density profile set by the progenitor optical depth, 
p = R(t v = \)/nr , and the X-ray opacities of Eqn. |2T] The escape fractions were computed for fixed wind 
densities and velocities of v. v =1100 and 560 km/s, respectively. Let fx.abs be the fraction of the X-rays 
energy which is absorbed. We also apply the correction from Eqn. [13] for the balance between cooling and 
expansion losses. Thus, the shock energy absorbed in the wind is 

Lx,abs = Lsfx,abs [l +texp/ t cool] ■ (A6) 

In the simple model, we just treat this in the same manner is the radiation from the optical transient, where 
it contributes 

LuAt) = L x . ahs B u (0, T(t))e- Tv ^ X) (A7) 



by direct emission and with the same fraction absorbed by dust. The shoc k model discussed in 34.41 uses a 



list of emission lines and the fraction of energy radiated in each line from lAllen et all (|2008) as the model 
for the spectrum produced by the absorbed X-rays, and computes the extinction and fraction absorbed based 
on the model line spectrum. Finally, the shock breakout luminosity was set by 

L P eak= 16nrlQ rat aTj est (A8) 

with Tci es t = 1500 K, heating the dust at radius ro to Tdest- We used a fixed temperature of T pea k = 50000 K 
and gave it a duration comparable to the light crossing time of the progenitor r*/c with r* = 5 x 10 13 cm. 
We only consider the "dust echo" from the breakout peak and ignore the short UV flash. 
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The dust luminosity is then 

L d (t) = (L peak (t)+Lo(t)+L x>abs (t))(l -f ,esc) (A9) 

and the dust temperature is 

Td= { l6na Qrat r d (t)i ) ' (A10) 
including the factor Q rat (average ratios of Planck factors) to make it easier to match this model to the 
DUSTY results. We then used a modified Planck function for the SED of the dust emission, 

L uM = L d B u (l,T d ). (All) 

We then include the effects of light travel time ("dust echoes") on the dust emissions, spreading the contri- 
bution of the emission from any time t uniformly over the interval t <t' <t + 2r(f) / c, so the total emission 
is 

L v =L v , {t)+L VtK {t) + {L v4 {t')). (A12) 

We also experimented with including direct emission from the progenitor, modeled as a black body with 
the luminosities from §|2] and a temperature of 3000 K. The models were then fit to the light curve data, 
with priors to match the parameters of the DUSTY models for R(tv = 1) in the progenitor models and the 
sequence of luminosity, optical depth and dust radius estima tes in the transient models. The shock velocity 



was constrained to match the estimates by lSmith et al.1 (120101) to 10% and the "Planck" term was constrained 
to have log Q rat = to 10%. 
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Fig. 1. — The spectral energy distributions of the SN 2008S (left) and NGC 300-OT (right) progenitors 
modeled as DUSTY winds. The underlying star is assumed to be a cold, = 2500 K, extreme AGB star, 
(badly) modeled with the black body SED shown by the dashed lines. Note, however, that no particular 
stellar temperature is required by the models because of the high optical depths. The solid curves show the 
DUSTY wind model for the observed SED assuming an inner edge dust temperature of 1500 K. The models 
have visual optical depths of ry = 300 and 750, re s pective ly, and stellar luminosities of L* = 10 4 6 L(p and 



lO^Lm. The SN 2008S data are fromlPrieto et all (12008b). the NGC 300-OT mid-IR data are from 



Prieto 



(2008a) (also Thomps on et all (120091)) and the NGC 300-OT optical limits are from iBerger & Soderberg 
J2008I) (also feerger et al.ld2009l) and bond et all J2009h ). 
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Fig. 2. — The stellar luminosity and wind density parameters for DUSTY models of the progenitors of 
SN 2008S (left, black) and the NGC 300-OT (right, red). These include all models roughly consistent 
with the photometry and having inner edge dust temperatures > 1000 K. The cyan curve shows the wind 
density predicted from the radius R(jv =1) where the visual optical depth is unity and a dust optical depth 
of ry = 137 cm 2 /g. This is for DUSTY's default gas-to-dust ratio of r gi / = 200, and the wind density can be 
raised as M/v w oc r g d by reducing the dust opacity per unit mass, Ky oc r"j, or vice versa. 
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Fig. 3. — The radii R(ry = 1) and R(tj = 1) where the visual optical (upper points) and Ifim (lower points) 
optical depth is unity for the SN 2008S (left, black) and the NGC 300-OT (right, red) progenitor winds. The 
vertical lines show the model likelihood- weighted mean luminosities. The cyan line through the R{ry = 1) 
points shows the simple fits given in Eqn.|2j while the line just below the R(tj = 1) points show the prediction 
for that radius given a simple p oc l/r 2 wind model. 
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Fig. 4. — The evolving SEDs of SN 2008S. Filled squares show measurements at that epoch. Filled triangles 
show upper bounds either from that epoch or a plausible bound set by the other epochs. Open triangles show 
plausible lower bounds set by the other epochs. The open pentagons show the SED of the progenitor star 
excluding the optical upper limits for clarity. The solid line shows the probability-weighted mean SED of 
the DUSTY models and th e dashed line shows th e dispersion of these SEDs. Epochs are relative to the same 
reference date as used by lBotticella et al.l (I2009T) . 
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Fig. 5. — The evolving SEDs of the NGC 300-OT. Filled squares show measurements at that epoch. Filled 
triangles show upper bounds either from that epoch or a plausible bound set by the other epochs. Open 
triangles show plausible lower bounds set by the other epochs. The open pentagons show the SED of the 
progenitor star excluding the optical upper limits for clarity. The solid line shows the probability-weighted 
mean SED of the DUST Y models and the dashed line shows the dispersion of these SEDs. The optical limit 
from lPrieto et al. ( 2010 ) is applied to the last two epochs but lies below the lower edge. Epochs are relative 
to 17 April 2008. 
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Fig. 6. — Simple light curve models for SN 2008S. The top, middle and bottom panels show the mid-IR 
(IRAC 3.6, 4.5, 5.8, 8.0 and MIPS 24/im), near-IR (J, H and K), and optical (UBVRI) light curves along 
with the simple model fits. Measurements are shown by filled squares, upper bounds by triangles. The 
late time optical/near-IR light curves primarily illustrate the time scales on which the absorption exterior to 
the shock becomes negligible - while correct for overall energetics, they are wrong in detail because they 
(incorrectly) assume that the SED of the shock emission is the same as the black body SED of the optical 
transient. 
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Fig. 7. — Simple light curve models for the NGC 300-OT. The top, middle and bottom panels show the 
mid-IR (IRAC 3.6, 4.5, 5.8, 8.0 and MIPS 24/zm), near-IR (J, H and K), and optical (UBVRI) light curves 
along with the simple model fits. Measurements are shown by filled squares, upper bounds by triangles. 
Models allowing modest changes in the transient temperature better fit the "kink" in the optical light curves. 
The same caveats apply here as in Fig. [6] The poorer fit for the R band is probably due to Ha emission. 
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Fig. 8.— The luminosity evolution of SN 2008S (left) and the NGC 300-OT (right). The points show the 
luminosity estimates from the DUSTY models. The solid curve shows a simple exponential plus constant 
model for the luminosity, while the heavy dashed curve shows the model from the light curve fits. The 
dashed curve shows the luminosity of the progenitor. 
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Fig. 9.— The optical depth evolution of SN 2008S (left) and the NGC 300-OT (right). The points show the 
ry estimates from the DUSTY models, while the heavy dashed curve shows the results from the light curve 
fits. 
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Fig. 10.— The dust radius evolution of SN 2008S (left) and the NGC 300-OT (right). The points show 
the radius of the inner edge of the dust from the DUSTY models while the heavy dashed line shows the 
model from the light curve fits. The right axis shows the visual optical depth scale of the progenitor wind, 
where the horizontal dashed lines mark the radii where the visual (ry = 1) and 7/im (77 = 1) optical depths to 
infinity were unity. The solid lines show the small grain dust destruction radius Rdest from Eqn. [T4lpredicted 
from the luminosity model in Fig. [H and the shock radius R s hock = v s t assuming velocities of v s = 1000 
and 600 km/s for SN 2008S and the NGC 300-OT, respectively. The dotted line shows the light radius 
Rlight = ct. Epochs above this line will be affected by finite light crossing times, although not significantly 
on the logarithmic scales of Figs.l8ll9land[T0l 



35 




10 1 



10 1 



-500 



o@@oo 



0) 

o 

s 



50% surviving opacity 



'^ r n\)i^ ieb 



6^ 




1 shock 



500 

time 



1000 

days 



1500 



2000 



Fig. 11. — Physical scales relevant for dust formation as a function of time and scaled for SN 2008S. The 
heavy solid line shows the radius Rf orm (Eqn. [Pfl) interior to which small grains would be destroyed. The 
heavy dashed line shows the radius of the outgoing shock. The horizontal dashed line shows the radius 
at which 50% of the dust opacity would be destroyed given a shock break out luminosity of 4 x 10 10 Lq 
(Eqn. [5]). The dotted shows the radius interior to which hydrogen can have recombined. Finally, the light 
black lines show the radii interior to which grains can have regrown to 0.003 /im, 0.01/im, and 0.03 /xm, 
respectively. These would move outwards in radius by a factor of two for the coagulation model. The 
cartoons at the top roughly illustrate where dust is located and its evolution in radius and optical depth. 
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Fig. 12. — Simple models for the reformation of dust in the wind of SN 2008S (solid curves). The lower 
curve assumes the growth time scale of Eqn. Q3] and the upper curve assumes a rate 4 times faster, corre- 
sponding to growth by coagulation rather than monomer accretion. The points and dashed curve show the 
optical depth estimates from the fits to the SEDs and the lightcurves, respectively, as in Fig. [9] Given the 
simplicity of the models, only qualitative agreement is to be expected. 
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Fig. 13.— The predicted X-ray properties of SN 2008S (solid) and the NGC 300-OT (dashed). The 
bottom panel shows the estimated luminosity including both absorption and the crude estimate of the balance 
between radiation and expansion. The middle panel shows the mean X-ray energy after absorption. The 
top panel shows estimates of the count rates for SWIFT (lower count rates, black) and Chandra (higher 
count rates, red). The points show the count rates for finding 3 counts in archival SWIFT (black diamond 
for SN 2008S open triangles for the NGC 300-OT) and Chandra (filled triangle for the NGC 300-OT) 
observations given their integration times. Multiple SWIFT epochs were combined into a single limit. 
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Fig. 14. — Alternate light curve models for SN 2008S. As in Fig. [6l the top, middle and bottom panels 
show the mid-IR (IRAC 3.6, 4.5, 5.8, 8.0 and MIPS 24/im), near-IR (J, H and K), and optical (UBVRI) 
light curves. The solid lines show the light curves using an emission line model for the shock spectrum, and 
the dashed lines add an unobscured 3000 K black-body with the luminosity of the progenitor. As expected, 
the differences between the models appear in the near-IR and optical as the optical depth to the shock front 
diminishes. The shock model has relatively little near-IR emission, so the near-IR emission never recovers, 
while the shock plus progenitor model is brighter in the near-IR than the original model because of the 
cooler progenitor black body temperature. The latter model assumes that there is no significant dust optical 
depth interior to the shock. 



